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Abstract 

Base station cooperation improves the sum-rates that can be achieved in cellular systems. Conven- 
tional cooperation techniques require sharing large amounts of information over finite-capacity backhaul 
links and assume that base stations have full channel state information (CSI) of all the active users in 
the system. In this paper, a new limited feedback strategy is proposed for multicell beamforming where 
cooperation is restricted to sharing only the CSI of active users among base stations. The system setup 
considered is a linear array of cells based on the Wyner model. Each cell contains single-antenna users 
and multi-antenna base stations. Closed-form expressions for the beamforming vectors that approximately 
maximize the sum-rates in a multicell system are first presented, assuming full CSI at the transmitter. For 
the more practical case of a finite-bandwidth feedback link, CSI of the desired and interfering channels 
is quantized at the receiver before being fed back to the base station. An upper bound on the mean 
loss in sum rate due to random vector quantization is derived. A new feedback-bit allocation strategy, to 
partition the available bits between the desired and interfering channels, is developed to approximately 
minimize the mean loss in sum-rate due to quantization. The proposed feedback-bit partitioning algorithm 
is shown, using simulations, to yield sum-rates close to the those obtained using full CSI at base stations. 
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I. Introduction 

Base station cooperation is an effective strategy to increase data rates and reduce outages in multiple- 
input multiple-output (MIMO) cellular systems [1]— [16]. Cooperative encoding at the base stations 
can be used to combat co-channel interference (CCI), paving the way for more aggressive frequency 
reuse, which can lead to higher data rates and simpler network configurations. Consequently, cooperative 
transmission is being considered for upcoming cellular standards like 3GPP long term evolution advanced 
that are targeting universal frequency reuse [17]. Base station cooperation entails sharing control signals, 
user propagation channel information and/or precoding data via high-capacity wired backhaul links to 
coordinate transmissions [1]— [16]. Most of the literature on cooperative techniques assumes that base 
stations have full channel state information (CSI) of all active users in the system. This is infeasible as 
receivers use finite-bandwidth channels to feedback CSI to base stations (in frequency division duplex 
systems) [18]. Another common assumption made in multicell cooperative literature is that base stations 
are connected using ideal backhaul links. In practical systems, however, the prohibitive cost of connecting 
all the base stations in the network using high-capacity links restricts backhaul capabilities [1]. While 
sharing a larger amount of information among base stations can improve performance, the extent of 
cooperation among base stations is restricted by the capacity-limited backhaul. Hence, it is important to 
develop cooperative techniques with limited feedback that maximize performance gains while ensuring 
a manageable load on the finite-capacity backhaul. 

Perhaps the least aggressive form of cooperation is joint resource allocation or joint scheduling by base 
stations in adjacent cells [2], [3]. These strategies typically result in a small load on the backhaul link 
and possess comparatively low complexity. For example, in dynamic and fractional frequency reuse, base 
stations exchange control-level information to cooperatively assign different frequency bands of operation 
to users in adjacent cells [2]. By utilizing the available spectrum more efficiently, the techniques in [2] 
yield higher sum-rates as compared to static frequency reuse, which does not account for varying user 
traffic. Alternatively, joint inter-cell scheduling assigns different transmission cycles to users in adjacent 
cells using opportunistic scheduling [3], introducing an expanded multiuser diversity gain relative to static 
frequency reuse. Frequency reuse and inter-cell scheduling do not utilize all the available frequency and 
time resources, respectively and hence, do not realize the performance gains that can be potentially 
obtained using base station cooperation [3]. 

Full cooperation among base stations leads to the highest sum-rates, at the cost of increased overhead 
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associated with the exchange of a greater amount of information among base stations. The sum-capacity 
maximizing solution was derived and found to be achieved by the dirty paper coding (DPC) [4]-[6], 
[19]. DPC requires joint precoder design, where all the base stations have perfect knowledge of the 
interference seen by each user. DPC was proposed for single antenna cellular systems in [4]. Multicell 
DPC using MIMO was investigated in [5], and an upper bound was derived on the sum-rates achievable 
using base station cooperation. While the approaches in [4], [5] assumed a sum-power constraint across 
all the base stations, DPC with more realistic per-base power constraints was considered in [6]. Multicell 
DPC is difficult to implement in practice due to the requirement for the base stations to have a non- 
causal knowledge of the interference [5], [7]. As a result, several low-complexity and relatively more 
practical joint transmission strategies were proposed in [5], [8] using full cooperation and per-base station 
power constraints. These sub-optimal linear techniques, including zero forcing, minimum mean square 
error or null-space decomposition, yield sum-rates that are far from that of multicell DPC due to the 
transmit power inefficiency introduced by the precoding matrices, especially for low -rank channels [5], 
[8]. Though the sub-optimal methods reduce the encoding complexity, they still require full CSI at the 
base stations and involve the exchange of a large amount of information among base stations, resulting 
in a prohibitive load on the finite-capacity backhaul links [1], [5], [8]. 

Partial cooperative strategies, where base stations exchange only the CSI of active users, offer a fair 
balance between ensuring a reasonable load on the backhaul links and attaining the performance gains 
using cooperation [1]. The shared CSI can be used by the base stations to design individual precoding 
matrices (or beamforming vectors, for single-stream transmission) on site to transmit exclusively to users 
within their own cell [4]-[6]. Most of the existing literature [9]— [15] also assumes that users can feedback 
full CSI to the base stations. A distributed approach for beamforming was studied in [9], [10] where 
transmit symbol vectors were designed using linear minimum mean square error (MMSE) estimation 
techniques. The forward-backward and the sum-product Kalman smoothing algorithms were used to 
iteratively obtain the transmit vector as an optimal linear MMSE estimate for a linear and a 2-D array of 
cells, respectively [9], [10]. The convergence of these iterative approaches is, however, not guaranteed. 
A distributed zero-forcing beamforming strategy was proposed in [11] using a scheduling algorithm to 
select users with the best desired channels in each cell. The sub-optimal approach is not feasible for 
practical setups as it uses a sum-power constraint across all base stations, and only satisfies an equal 
power per base station asymptotically for a large number of users per cell [11]. Distributed transmit 
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beamforming was also investigated in [12], where the authors propose an iterative algorithm to minimize 
the transmit power that does not necessarily maximizie the sum-rates. The authors of [13]— [15] propose 
a beamforming solution to approximately maximize the sum-rates for only a two-cell system. 

In this paper, we propose a linear beamforming strategy to approximately maximize the sum-rates 
at high signal-to-interference noise ratios (SINR) in multicell systems, using partial cooperation. The 
proposed beamforming method is similar to the two-cell approaches in [13]— [15]. The cellular setup 
considered is based on the Wyner model [20], where neighboring base stations share only the CSI to 
beamform independently to one active user in the cell with a single (dominant) co-channel interferer. 
The proposed algorithm is also applicable to a finite linear array of cells. The variation in the strength 
of the desired and interfering signals is taken into account to maximize the sum-rate using a generalized 
eigenvector approach. We first assume full CSI availability at base stations and a high-capacity backhaul 
to arrive at the closed-form expressions for the mutlicell beamforming vectors. We then consider the more 
realistic limited feedback and capacity-limited backhaul in the rest of the paper. The beamforming strategy 
presented is non-iterative in nature and hence, does not have convergence issues unlike the solutions in 
[9], [10]. The linear approach in this paper uses explicit per-base station power constraints, in contrast 
to distributed zero-forcing beamforming in [11]. Further, it possesses low-complexity, and results in a 
smaller burden on the backhaul link, in comparison to the full cooperation strategies in [4]-[6], [8], while 
yielding sum-rates that are reasonably close to those of multicell DPC. Also, it can be applied to models 
with any number of linearly arranged cells, and not just two-cell situations as in [13]— [15] . 

Implementing the linear beamforming strategy requires full CSI at the base stations to design the 
transmit beamforming vectors, which may be impractical. Further, all the strategies described in [1], [4]- 
[16] require full CSI at the base stations. In practice, as feedback links have finite bandwidth, quantized 
CSI is fed back to the base stations using the concept of limited feedback [18]. While limited feedback 
for single-cell systems is well researched (refer to [18], [21] and the references within), comparatively 
less work has been done in the multicell scenario. In contrast to the single-cell case, the CSI of multiple 
channels needs to be fed back for the cooperative-based strategies. Hence, it is not even clear how to 
apply the existing multiuser limited CSI feedback techniques [22], [23] to multicell systems. Also, the 
varying strengths of the desired and interfering channels need to be taken into account when developing 
a multicell limited feedback strategy that will efficiently utilize feedback resources. Single-cell multiuser 
limited feedback does not have to deal with the variability in signal strengths, due to the single-channel 
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feedback. Note that the existing work in multicell cooperation that considers the difference in signal 
strengths of the desired and interfering channels [8], [11], [16] assumes full CSI at base stations. 

In this paper, we develop a limited feedback strategy for the proposed beamforming algorithm that will 
partition bits between the desired and interfering channels as a function of their relative strengths. For 
analytical reasons, the desired and interfering channels are quantized using random vector quantization 
(RVQ), i.e. the quantization vectors are independently chosen from the isotropic distribution on a unit 
hypersphere [24], [25]. We first present a simple model for limited feedback of CSI in multicell systems. 
The model requires each user to feedback quantized CSI to its own base station. Adjacent base stations 
are connected by a backhaul link, which is used to transfer the quantized CSI of only the interfering 
channels. Hence, the requirement of a global high-capacity backhaul link is eliminated. We quantify the 
reduction in sum-rate caused by quantizing CSI using RVQ by deriving an upper-bound on the mean 
loss in sum-rate. We also present a feedback-bit allocation technique that approximately minimizes the 
mean loss in sum-rate and hence, utilizes the available feedback resources efficiently. Using the proposed 
algorithm, we show that the sum-rates approach those of multicell DPC even for a limited feedback 
scenario. 

We first extend the two-cell cooperative beamforming approach in [13]— [15] to a multicell system where 
adjacent base stations exchange only CSI to approximately maximize sum-rate. The main contributions 
of this paper are as follows. 

• We propose a novel and simple limited feedback model for the multicell system, where each user 
feeds back quantized CSI of the desired and interfering channels to its own base station. In the 
proposed model, the neighboring base stations exchange only the quantized interfering CSI. This 
makes the load on the backhaul link manageable. By ensuring that only the adjacent base stations 
are connected, the requirement of a globally connected backhaul is eliminated. 

• We analyze the performance of the beamforming technique presented in this paper with limited 
feedback by deriving an upper bound for the mean loss in sum-rate due to CSI quantization using 
RVQ. It is shown via simulations, that the upper bound is reasonably tight. 

• We present a new feedback-bit partitioning algorithm to allocate the available feedback bits between 
the desired and interfering channels to approximately minimize the mean loss in sum-rate due to 
limited feedback. 

This paper is organized as follows. We first describe the system models used in this paper in Section 
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PH We then present the multicell beamforming strategy, and compare it with two existing beamforming 



techniques for the full CSI case in Section III In Section [TVl we consider the limited feedback scenario 



for the beamforming algorithm in Section III We propose a feedback-bit allocation strategy to partition 
the available bits between the desired and interference channels in Section [V] The simulation results 
verifying the accuracy of the limited feedback algorithm are presented in Section [VT] The concluding 



remarks are provided in Section VII and detailed proofs are presented in Appendices [I] - [IV] . 



Notation: In this paper, X refers to a matrix and y stands for a vector. The transpose and conjugate of 
X are given by X T and X c , respectively. The Hermitian transpose of X is given by X*. The inverse 
and pseudo-inverse of X are given by X -1 and X^, respectively. An identity matrix of size R x R is 
denoted by I/j. refers to the expectation. ||x|| stands for the Frobenius norm of x. J\f c {^,cr) refers 
to a complex Gaussian random distribution with mean \i and variance a. 

II. System Model 

Consider the multicell setup shown in Fig. [T] for an array of linearly arranged cells. The multicell 
system used in this paper is based on the Wyner model [20]. We assume that the base station in each 
cell serves a single active user [4], using intra-cell time division multiple access. Each user is assumed 
to face interference from one of its neighboring cells, as shown in Fig. [T] The received signal power 
of the desired and interfering signals is a function of the user's location in the cell. A similar approach 
was adopted in [8], where a user at the cell center receives a signal from only its corresponding base 
station^] while a user at the cell edge is subjected to interference from an adjacent cell as well. The 
results in this paper can also be applied to modifications of the Wyner model, including the circular 
arrangement of cells proposed in [11], and a finite linear cellular array. Note that the circular extension is 
a generalization of the two-cell case and is equivalent to the Wyner model as the number of cells becomes 
large. The generalized eigenvector beamforming strategy proposed in this paper can be extended to the 



finite-array case by accounting for the edge effects, as will be described in Section III The Wyner model 
and its extensions are widely used in literature as they model the central factors of a cellular system like 
fading and inter-cell interference, while retaining analytical tractability [8], [11], [16]. We assume that 
the backhaul link used for information exchange between base stations is ideal and that the time delay 
associated with feedback and cooperation is zero. 

'in our paper, we assume that the interference from the adjacent cell is negligible at the cell center. 
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The number of cells in the multicell system is denoted by K, where K goes to infinity for the Wyner 
model. We index the users in each cell by the base station they obtain their desired signal from, i.e. the 
k th base station services the k th user, for k = 1 , . . . , K. We assume that all the base stations are equipped 
with Nt antennas, while each user supports a single receive antenna (i.e. multiple-input single-output, or 
MISO, system). The channel corresponding to the desired signal between the k th base station and k th 
user is denoted by € C tXl . The interfering channel between the k th user and the {k + l) th base 
station is given by g^+i G C NtXl . This is illustrated in Fig. [I] The symbol transmitted from the k th base 
station (intended for the k th user) is denoted by Sk, where the transmit power, E{|sfc| 2 } is normalized to 
one. The transmitted signals are subjected to large-scale fading, which includes distance-dependent path- 
loss and shadowing effects, and small-scale fading. After averaging over the small-scale fading effects, 
the desired and interfering signal powers received at the user terminal are denoted by 7fc,(d) an d 7jfc,(i)» 
respectively, for the k th user. To facilitate analysis in different loss settings, we let Jk,(i) = a klk,(d)' 
where a* £ [0, 1] (i.e. the interfering signal strength can at most be equal to that of the desired signal). 
Note that a similar parameter is used in [8] to model the SNR of the interfering signal with respect to the 
received signal. Using the narrowband flat-fading model, the baseband discrete-time input-output relation 
for the user in the k th cell is given b)H 



where yk € C is the received signal at the k user and G c w * is the beamforming vector at the 
k th base station. Finally, G C is complex additive zero-mean white Gaussian noise at the receive 
antennas, with E{|n/c| 2 } = N . 

The signal to interference noise ratio (SINR) of the k th user is given by 



Vk = ^/j k ,(d)K { kSk + y/TA^gfc+iffc+iSfc+i + n k , 



(1) 



SINR fc 



7fc,(d)|hfcffc 



(2) 




(3) 



afc|gfc+iffc+i| 2 + ^ 



Note that Pfc,(d) = 7 ^- (d) is the received SNR of the desired signal, which is independent of the beam- 
forming vectors. The base stations are assumed to have perfect knowledge of Pk t (d)- This is a popular 



2 We drop the discrete-time index for sake of convenience. 



assumption in literature [26]-[28]. It was shown in [28] that SNR quantization does not effect the sum- 
rates of a single-cell multiuser MIMO system signicantly. To the best of our knowledge, the effect of 
SNR quantization on the sum-rates of a multicell system has not yet been investigated. We analyze the 
impact of SNR quantization in our future work on cooperative beamforming. In this paper, however, we 
assume that the base stations have perfect knowledge of the SNR so we can concentrate on quantizing 
the channel direction. The sum-rate of all the users within the system, R s , is expressed as 

E s = ^log 2 (l + SINR fc ) . (4) 

k 

It is evident that the design of the beamforming vectors {fk\k=i influences the SINRs obtained at each 
of the users and hence, the sum-rate of the multicell system. The dependence of SINR^ on both, and 
ffc+i (as given in ([3])) implies that a joint optimization across all the active users in the multicell system 
is required to maximize the sum-rate. In this paper, we use a high SINR approximation to remove this 
interdependency of users and thereby avoid an explicit joint maximization. 

III. Designing Beamforming Vectors Assuming Full CSI 

In this section we present a multicell cooperative beamforming approach. First, we describe two popular 
existing beamforming strategies, non-cooperative eigen-beamforming and cooperative zero-forcing beam- 



forming, which we use for comparison with the proposed algorithm in Section VI Eigen-beamforming 
maximizes the desired signal strength, while zero-forcing beamforming nulls out the interference. Next, we 
describe the proposed generalized eigenvector based beamforming strategy that approximately maximizes 
the sum-rate at high SINR, in a multicell system. We consider the multicell setup given in Section |TT] and 
assume that full CSI is available at the base stations. 

A. Eigen-Beamforming and Zero-Forcing Beamforming 

In eigen-beamforming, each base station transmits in a non-cooperative manner by ignoring the co- 
channel interference from neighboring base stations, as in the case of single-cell single-user beamforming 
[29]. For MIMO systems, the beamforming vector is chosen to be the eigenvector corresponding to the 
maximum eigenvalue of the MIMO channel. For a MISO system, the eigen-beamforming vector, f^, is 
reduced to 

h = A ■ (5) 



Note that (|5]) maximizes the numerator in ([3]). Eigen-beamforming maximizes the data rate that can be 
obtained in single-cell single-user systems, with no inter-user or inter-cell interference. 

The zero-forcing approach in [5] requires base stations to share data in addition to CSI. To ensure 
fair comparison with the beamforming strategy presented, we implement zero-forcing with single-cell 
processing in a way that only CSI is shared between base stations, using a principle similar to that in [5]. 
Adjacent base stations exchange the interfering channel state information between each other, so that the 
k th base station has the knowledge of both and g*.. Assuming full CSI at the transmitter, the multicell 
zero-forcing criterion requires that |g^ffc| 2 = 0, where fff*. = 1. Denoting = [hT, gTy, the zero 
forcing beamforming vector is designed by setting 

Note that g^Wfe(:, 1) = 0, and hence, the solution in ([6]) minimizes the denominator in ((3]). 

B. Proposed Beamforming Strategy for Approximately Maximizing Sum-Rate at High SINR 

In this section, we present a linear multicell cooperative beamforming strategy (for the setup in 
Section [II]), based on the two-cell beamforming solution proposed in [13]— [15]. The approach in [13], 
[14] approximately maximizes the sum-rate for a two-cell MISO system where base stations are assumed 
to have full CSI. In this paper, we extend the analysis to develop a cooperative beamforming technique 
that approximately maximizes the sum-rates in the high SINR regime in a multicell system. 

At high SINR, since log(l + SINR) » log(SINR), Q can be expressed as 

fl a «^log 2 (SINR fc )=log 2 (jJSINR Jfc ) . (7) 

k \ k ) 

Thus, maximizing the sum-rate, R s , approximately at high SINR involves maximizing the product of the 
SINR^] The beamforming vectors {fk\k=i are found by solving 



(k Ih^f P 
TT , T ' k , 2 3-1 s.t. ||f k || 2 = l, for k=l,...,K . (8) 

Taking advantage of the commutativity of the multiplication operation, the objective function can be 



3 We show, using simulations, in Section 



VI 



that the approximation is tight for small au or large p^. 



(d)- 
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written as 

/ |h T ft,| 2 \ 
R s>opt = max VlogJ , ' fc , 2 — s.t. ||f k || 2 = 1, k = l,...,K. (9) 

The A; th base station uses CSI of and to obtain the optimal if,, implying that the sum-rate can be 
approximately maximized by each base station exchanging CSI only with its neighbors. This eliminates 
the need for global CSI knowledge and reduces the load on the finite-capacity backhaul link. 

Now, ([9]) deals with maximizing a function consisting of K independent variables, {ffcjfLi- Hence, 
it can be split into K independent problems. Thus, the high SINR approximation is used to remove the 
interdependency of the users in Q. The optimal linear beamforming vector, fj. )0p t, is the solution to 

f*Rh fc f 1A 

iA-ont=arg max — — — -, (10) 

* P f:||f|| 2 = l f*Rg fc f 

where R hfc = Yi k h* k and R gfc = p fc _i,(i)gfcgfc + IjV t , and p k -i,(i) = afc-iPfc-l,(d) is the received SNR 
of the intereference signal at the (k — l) th user. The expression in (K)i is the well known generalized 



Rayleigh quotient [30]. Since Rh fc and R gfc are Hermitian and R gfc is positive definite^} the solution to 
( fT0| ) is given by the generalized eigenvalue decomposition: 

R hfc f = A fc R gfc f, (11) 

where denotes the eigenvalues of Rg^R^ [31]. Due to the rank deficient nature of Rh fc (= hfehp, 



there exists only one non-zero eigenvalue. The solution to ( |T0| ) will then be equal to the generalized 
eigenvector corresponding to the non-zero (maximum) eigenvalue. The solution for ffc i0pt is invariant to 
an angular rotation, 9, of the beamforming vector, e 3 ' ffc )0 pt- Hence, there can be infinitely many solutions 



to ( pLOp - The unit-norm and rotationally invariant nature of ffc )0p t implies that it is a point on the Grassmann 
manifold, used popularly in single-cell limited-feedback literature [18]. 

Note that while ( [TT] ) can be used to approximately maximize the sum-rate for both an infinite linear 
array (in Fig. [T]) and a circular arrangement of cells, it can also be used for a finite cellular array by 
accounting for edge effects and using ([TT]) for k = 2, . . . , K. The base station in the first cell of the array 
does not cause interference (ao = 0), implying that optimal beamforming vector is given by ([5]). Hence, 



(111 can be used for any K, while the solution in [13], [14] approximately maximizes sum-rates for only 



4 R gfc is positive definite (and hence, full rank) as x*R gfc x = x*(pfe„i i (i)gfcgfe + Ijv t )x > 0, for all 
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K = 2. Also, while it is recognized that the ratio in ([9]) is similar to the signal-to-leakage-noise -ratio 
(SLNR) proposed in [32] for multiuser MIMO, the single-cell SLNR solution does not bear a direct 
impact on the sum-rates obtained and is shown to yield low sum-rates in the high SINR regime [33]. 
In contrast, we show in ((9]) that (approximately) maximizing the sum-rate at high SINR is equivalent to 
maximizing SLNR at each base station for the multicell setup in Section [II] 

IV. Designing Beamforming Vectors Using Limited Feedback 



In Section III we presented a beamforming approach to approximately maximize the sum-rate at high 
SINR in a multicell setting by assuming that full CSI of the desired and interfering channels was available 
at the base station. A reasonable way to inform the base station about the channel state is through a 
limited feedback channel. In this section, we propose a limited feedback approach for multicell generalized 
eigenvector beamforming, where the quantized CSI of both the desired and interfering channels is fed 
back over the feedback link. The available feedback resources can be utilized efficiently in the multicell 
case by taking into account the relative strengths of the desired and interfering channels. Note that this 
is in contrast to single-cell multiuser limited feedback case, where only the CSI of a single channel is 
fed back [22], [23]. 

The limited feedback model for the multicell setup in Section [TT] is described in Fig. [2] It is assumed 
that the k th user can perfectly estimate the desired and interfering channels, h k and g^+i, using separate 
training symbols from the k th and (k + l) th base stations, respectively. The k th receiver uses -Bfc,(d) 
and fljfe m bits to feed back quantized versions of the and gfe+i, respectively, over a limited feedback 
channel where B k ^ + B k n\ = 2? tot is fixed. Note that B k ^ and B k m will depend on the relative 
strength of the interfering and desired signals (a k ). For example, when a k — > at the cell center, the 
contribution of the desired channel towards the SINR is greater than that of the interfering channel and 
it becomes more important to reduce the loss due to quantization of as compared to that due to gfc+i- 
This can be done by allocating most of the feedback bits to quantize h^, i.e. Sjt (a) -Btot- Hence, the 
two channels, and g k +i, are quantized separately using two separate codebooks of (variable) sizes 
2 B M d > (denoted by W fc (d ))) and 2 Bk ^ (given by W k ^), respectively. 

The unit-norm direction of the estimated desired (and interfering) channel is given by = hfc/||hfc|| 
(gfc+i = gfc+i/||gfe+i||)- The channel directions are then mapped to the respective codebook entries with 
the smallest angular separations. The quantized vectors for the desired and interfering channels, denoted 
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by hjt and gk+i, respectively, are obtained using 

hfc = arg max |h£w| 2 = arg max cos 2 (Z(hfc,w) ) , and (12) 

wew fc ,(d) wew fc ,( d ) \ / 

gfc+i = arg max |g^. +1 w| 2 = arg max cos 2 (Z(g fe+1 , w)) . (13) 

w6Wk,(i) w6W*,(ij 

Since the channel gains, ||h&|| and ||gfc+i||, are scalars and can be quantized easily, we concentrate on 
quantizing the channel directions, and gfe + i and assume that the base stations have perfect knowledge 
of ||hfc||. This is a popular assumption made in multiuser MIMO literature [26]-[28]. As quantizing 
channel gain is (to the best of our knowledge) not investigated for multicell systems, we use the perfect 
||hfe|| and ||gfc+i|| knowledge at base station assumption. We will analyze the impact of quantizing channel 
gain in our future work on limited feedback for multicell systems. 

From (|9]), it is seen that the k th base station requires information about both hfc and g& to compute f^. 
In the limited feedback model presented in this paper, the k th base station has information of both 
and g/t+i (in addition to perfect knowledge of ||hfc|| and ||gfc+i||)- The interfering channel information, 
||gfe+i||, is then sent from the k th to the (k + l) th base station via the backhaul link. In this manner, 
each base station has knowledge of not only its own desired channel, but also of the interference that is 
causing to the user in the adjacent cell. 

As codebook design for multicell systems is a topic of ongoing research, we use RVQ for channel 
quantization to facilitate analysis. If we use B bits for feedback, then each of the 2 B codebook vectors 
is independently chosen from the isotropic distribution on the N t dimensional unit sphere [24], [25]. 
The beamforming vector at the k th base station using limited CSI feedback, f& is then computed as the 
generalized eigenvector satisfying 

R^ffc = A fe Rg fc 4, (14) 

where Rj^ = [|hjfc|| 2 hfehj£ and Rg fc = Pfe-i^qllgfepgfcg^ +Ijv t . Quantization of CSI leads to a loss in 
the sum-rate. In this paper, we develop a feedback-bit partitioning strategy to approximately minimize 
the mean loss in sum-rate given by 

E{AR S } &E{R sMl - R S , LF } 
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where R Si m\ is the sum-rate obtained using the high SINR approximation with full CSI at the base station, 
as given in Q and R s ,lf refers to the sum-rate using limited feedback. Also, ffc+i,f u ii corresponds to the 



full CSI beamforming vector (from (111). 



V. Optimizing Feedback Bits to Minimize the Mean Loss in Sum-Rate at High SINR 



The beamforming strategy proposed in Sections III-B and IV requires channel state information of both 



the desired and interfering channel at the base station. As the number of available feedback bits, B to t is 
fixed, it is possible that allocating feedback bits between the desired and interfering channels can further 
improve the limited feedback sum-rate. In this section, we propose a feedback-bit allocation strategy 
to numerically evaluate the number of bits required to quantize the desired and interfering channels 
using RVQ. We model all the channels by the Rayleigh fading model, where each entry is a zero-mean 
unit-variance complex Gaussian independent and identically distributed (i.i.d.) random variable according 
to M C (0,1). While it is recognized that Rayleigh fading does not model realistic propagation channels 
accurately, we use the i.i.d. assumption in the limited feedback analysis to obtain closed-form expressions 
of the mean loss in sum-rate and feedback-bit partitioning algorithms. 

Maximizing the mean sum-rate using generalized eigenvector beamforming with limited feedback is 



equivalent to minimizing E{Ai? s } in ( p3] >, which is rewritten as 

E{AR s} = £ E (log 2 ( M^L 2 ) 1 + £E (log 2 ( Pk fff +l]2 + ' ) ) , (16) 
V 1 V l^ffcl 2 Jj k { \Pk,(i)W k+ ih+iMl\ 2 + 1 ) ) 



7V (d) T fc , ( i) 

by interchanging the terms inside the logarithms. The mean loss in sum-rate can be viewed to be 
a contribution of two terms, namely T k r d \ and T k n\, corresponding to the mean loss resulting from 
quantizing the desired and interfering channels, respectively. Note that T fc ( d ) and T k depend only on 
-Bfc,(d) an d Bk,(i)> tne num ber of bits used to quantize the desired and intereference channels, respectively, 
and are independent of TJ ( d ) and T^m, for k ^ I. Therefore, minimizing ( [To] ) is equivalent to minimizing 
Afe = T k _( d ) + T fc (i) for each k. Since obtaining a closed form expression for ( [To] ) is complicated, we 
derive an upper bound for E{Ai? s } in terms of bounds on TW d ) and T k ^y 



Before we proceed, we present a lemma that is used to derive an upper bound on (|T6j). Let H denote 
the set of all unit-norm channel directions, h. 
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2 



Lemma 1: The mean of log 2 (^cos |^Z(h,h)JJ (in (12l) is given by 

E W ,w{iog 2 (co S 2 (z(h,h)))} = io g2 ( e )^ (-ir E ?■ 

i=0 ^ ' Z=l 

Proof: The proof is given in Appendix [EJ ■ 
The following theorems present upper bounds on Tfcfd) an d Ifc (n. 

Theorem 2: The mean loss in sum-rate due to quantizing desired channel of the k th user, T k in 



( [16] ) is upper-bounded by 



r fej(d )<-iog 2 ( e ) e k-ir £ 

_L=o v 7 z=i 
Proof: The proof is given in Appendix [II] ■ 

Theorem 3: The mean loss in sum-rate due to quantizing interfering channel of the k th user, T k ^ in 



( fT6| ) is upper-bounded by 



<log 2 (l + p fei(i) iV t 2^ «/3^2 



7 Nt-1 



Proof: The proof is given in Appendix III 



Using Theorems [2] and |3j the difference in the mean sum-rates assuming full-CSI and using limited 
feedback, E{A7? S }, is upper bounded by 

E{AR S } < E^^ + Pfc.O)^ 2 ^ 1 ^^^^^)) 

- Elofe(e) E ( \ )(-!)' E }■ <") 

fc i=0 V y Z=l 



For a system with two transmit antennas at each base station, we substitute N t = 2 in (17 1 to get 

E{Ai? s }<Elog2(l + 2^, (i )2^"»/3(2 B M, 52 ))_^i g 2(e) ^ / " (-irEy 
fc fc i=0 ^ 2 ' 1=1 

Using series expansions, it can be shown that E^ 2 = o' (d) C ( — 1)* 2_w=i 7 = — 2 _Bfc (d) . Further, 
/5 (2 Bfc -< i >,2) = l/(2 B *.w(2 Bfc .« + 1)) and B fci(i) = 5 to t - #fc,(d)- Hence, (p) is rewritten as 

E{AR S } < Elog 2 (l + 2p m 2Btot J k {d) - 1 ) + 2'^ log 2 (e) . (19) 

fc s ^ , 

To find the optimum number of bits to quantize the desired and interfering channels, we use the relation 
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B k (;) = i? to t — -^fc,(d) i n '19). For simplicity, we denote the right hand side of expression ( fT9"] ) by A k . 
Treating B k ^ as a real number, we first show in Theorem |4] that A^ is convex in B k ^ £ [0,i?tot] ; 
and then use the result to compute the optimum number of desired and interfering feedback bits. 



Theorem 4: The minimum value of the upper bound, A k in ( [19] ), is obtained at the value of B k ^ 
equal to 

£fc,(d) real = Stat " log 2 (l + P k ,(i) + V W,(i)2 B - +1 + (p fc , (i) ) 2 ) , 

where -Bfc,(d) real S [0, B tot ] is a real number. 

Proof: The proof is given in Appendix [IVJ ■ 

Note that the minimization in Theorem [4] is over the set of real values and, hence, B k ^ real is not 
necessarily an interger. Since the upper bound A k is convex in -B fej ( d ) rcal , we only need to consider the 
ceiling and floor of -Bfc,(d) real (denoted by LSfc,(d) real J anc ^ ( d ) real ~|, respectively) to find the optimal 
number of bits to quantize the desired channel [34], [35], as given in Corollary [5] 

Corollary 5: The optimum number of desired and interfering feedback bits are given by B k ,(d)° pt an< i 
B k , (i) opt = Btot - £ fc ,(d)° pt respectively, where B k ^ is either [B kM) ^ l \ or [i? fe , (d) real l . 

Since Pk,(i) = Pk,(d) a k> it is clear from Theorem |4] that B k ,(d) real ( an d hence, -Bfc,(d)° pt ) i s a function 
of Pfc^d), ot k , and -Btot- Note that for a k = 0, B k ^ opt = B tot . This implies that when the interfering 
signal has zero power (or there is no interferer transmitting), all the available bits will be used to quantize 
the desired signal. In the absence of interference, it is evident that the strength of the received signal, 
Pk,(d) does not affect the bits allocated. In contrast, when a k = 1, B k ,(d)° pt increases as p k ^d) decreases. 
This makes intuitive sense because if the desired signal strength is high, fewer feedback bits are assigned 



to the desired channel and vice versa. In Section VI we show that the values of (B k ^,B k ^) from 
Theorem 0] match the simulation results. 

VI. Simulation Results 

In this section, we present simulation results for the case when full CSI is available at the base station 
and then present results for the limited feedback scenario. We show that the feedback-bit allocation 
strategy described in Section [V] approximately minimizes the mean loss in sum-rate due to quantization. 
Simulations are also used to verify that the partitioning of feedback bits between the desired and interfering 
channels proposed in Theorem [4] and Corollary [5] matches numerical results. 

Unless otherwise stated, we assume that all K users have the same received desired and interfering 
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signal strengths, i.e., Pk,(d) = P(d) ar *d a k = ol, for all k, for simplicity. We provide simulation results 
for the asymmetric case at the end of this section. 

A. Full CSI Case 



In Section III a high SINR approximation was made to approximately maximize the sum-rates in 
a distributed fashion across base stations. In Fig. [3] we verify that the approximation yields sum-rates 
that are reasonably close to the actual values for p^ values as small as 5 dB, for Nt = K = 4 and 
a = {0.001,0.1,1}. In Fig. |4j we show that as K (number of cells or users) increases, GEBF yields 
sum-rates that are very close to the upper-bound using multicell DPC for K users for N t = 4 and 
/9(d) = 10 dB. Though the difference in sum-rates is maximum for a = 1 because of a general drop 
in the SINR of each user due to increased interference from adjacent cells, it is seen that the multicell 
GEBF presented in this paper yields sum-rates reasonably close to the upper bound from DPC. 

B. Limited CSI Feedback 

The upper bound on the reduction in sum-rate (due to quantization) was derived in Section |Vj by 
assuming that the number of bits available for both, the desired and interfering channels are large enough 
to ignore quantization errors. Hence, in the simulation results presented in this section, for a B tot bit 
feedback system, we assume that S/t (a) £ [3, -Btot — 3], as using RVQ codebooks of size 2 3 = 8 
leads to low values of quantization error. In Fig. |5J for p( d ) = 10 dB, B to t = 15, K = Nt = 2 and 
a = {0.001, 0.1, 1}, it is shown that the upper bound derived in ([19]) is reasonably tight near the minimum 
mean loss in sum-rate. It is evident that partitioning adaptively the feedback bits, as a function of the 
location of the mobile terminal, is an effective method to increase the sum-rates obtained in limited 



feedback systems. This implies that the upper bound derived in (19 1 can be used to evaluate the number 



of quantization bits needed to minimize the mean sum-rate loss. Further, Fig. [5] clearly demonstrates the 



convexity of the actual loss in sum-rate and the upper bound in (jT9j). 

Next, we present simulation results to show that GEBF with the proposed feedback-bit partition- 
ing algorithm outperforms non-cooperative eigen-beamforming and multicell zero-forcing beamforming 



strategies with limited CSI feedback. Since the eigen-beamforming approach described in Section IV 
is non-cooperative, all the available feedback bits are used to quantize the desired channel only, i.e. 
-E>fc,(d) EBF = B tot , whereas for multicell zero-forcing approach, B k ^ ZF = B k ^ ZF = B tot /2 is 
assumed. Fig. [6] shows the sum-rates that can be obtained using the three transmission strategies for 
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a two-cell (two-user) scenario with B tot = 6 and pr d \ = 10 dB and K = N t = 2. It is seen that the 
proposed limited feedback algorithm yields higher sum-rates. As a increases, the sum-rates using both 
GEBF and non-cooperative eigen-beamforming reduce, but that of eigen-beamforming falls off more 
drastically. This is because as a — ► 1, interference becomes significant and the beamforming vector 
needs to reduce the interference in addition to increasing the desired signal strength. Zero forcing with 
full CSI yields a constant sum-rate for all values of a due to the nulling out of the interference term in 
the denominator of Q. 

The sum-rates obtained using the proposed bit-partitioning technique as B tot increases are compared 
in Fig. [7] for N t = K = 2 and p (d) = 10 dB for a = {0.001,0.1,1}. It is seen that the cell-edge 
users (a — » 1) requires a larger B to t as compared to users in the cell-center (a — > 1). This is due to 
the requirement of cell-edge users to quantize the interfering channel (in addition to the desired channel) 
with a sufficiently high resolution. In contrast, users in the cell-center that have weak interfering signals 
have to only ensure that the desired channel is allocated sufficient number of feedback bits. Hence, the 
feedback-bits, Btot can be varied adaptively as a function of the user location within the cell to yield a 
given sum-rate. 

In Fig. [8] we show the variation in the number of bits allocated to the desired and interfering channels 
as a function of the interfering to desired signal ratio, given by a = 101og 10 (a) dB. In Fig. [8j for 
K = N t = 2 and = 10 dB, we can see that when the difference in path-loss is greater than 37 dB, 
{B k ( d ), B k m) = (8,0). As the path-loss difference decreases, -Bfc,(d) reduces and B k ^ increases, due to 
an increasing need to quantize the interfering channel with greater resolution. When both the desired and 
interfering channels have the same signal strength, it is seen that (B k ^,B k ^) = (2,6), i.e. reducing 
the interfering signal strength is 'more important' to maximize the sum-rate than increasing the desired 
signal strength. 

Finally, we consider the case where the users are randomly located in their cells and have different 
interfering to desired signal ratios, i.e. a k ^ cti,k ^ I. Hence, the users can have different feedback 
bit allocations. The performance of the proposed feedback-bit allocation strategy with GEBF can be 
examined for the asymmetric scenario by evaluating the average data-rate in a cell, where the average 
is taken over all possible user locations in the cell. Assuming B tot = 6 at each user, we consider 1000 
cells (users) arranged in a linear fashion, with a k uniformly distributed between [—40 dB,0 dB], for all 
k. Adjacent cells cooperate by sharing the interfering channel CSI to compute the bit partitions for each 
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user independently. In Fig. |9j it is seen that the average data-rate in a cell using the proposed feedback-bit 
allocation strategy is reasonably close to that obtained using GEBF with full CSI. We also show that 
the data rate obtained using equal bit partitioning in GEBF, i.e. -Bfc,(d) = ^k,(i) = -Btot/2 for all k is 
less than that obtained using the proposed bit allocation, especially for larger p( d y This makes intuitive 
sense as the equal bit allocation does not reduce the quantization error, which grows with larger p( d ), 
while the proposed strategy allocates bit adaptively to minimize the overall impact of quantization error. 
We also plot the average data rates obtained using the non-cooperative eigenbeamforming and multicell 
zero-forcing approaches (described in Section III-A[ ) with limited feedback. The proposed feedback-bit 



strategy with GEBF outperforms both, EBF and ZF, especially at larger values of p^ d y Hence, it can be 
concluded that the proposed limited feedback strategy is effective in multicell systems. 

VII. Conclusion 

We first presented a beamforming strategy that uses cooperation among base stations in a multicell 
MISO system to approximately maximize the sum-rate at high SINR. A linear array of cells based on 
the Wyner model is considered, in addition to a circular extension of the model. For the limited feedback 
scenario, we presented a feedback-bit allocation strategy to approximately minimize the mean loss in 
sum-rate caused due to quantization using RVQ. The proposed technique relies on the relative strength 
between the desired and interfering signals within the cell to allocate bits adaptively to each of the two 
channels. A closed form expression for the number of feedback bits was derived. Using simulations, it 
was shown that the multicell beamforming approach yields sum-rates reasonably close to those obtained 
using multicell DPC, with full CSI. The average data rate in a cell using the feedback-bit allocation 
strategy was also shown to reasonably close to the full CSI case, verifying that the proposed multicell 
limited feedback algorithm yields high sum-rates using partial cooperation. 

Appendix I 
Proof of LemmaQ] 

It was shown in [36] that the probability density function of v = cos 2 (^(h, h)) is given by 

N i(N t -l) 

£ £ 

i=0 j=l 



N i(N t -l) 

«->=£ £ CD( ( 7 ^c-iw. 
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where N = 2 Bk - {i) is the size of the RVQ codebook. The mean, E{log 2 (^)}, is then computed as 



which is written as 



E{log 2 (z,)} = bg 2 (e) / hi{v)f v {v)dv, 
Jo 



i=0 ^ 1 ' .7=1 V J / J 

i(N t - m 



(20) 



E{log 2 (z/)} 



log 2 (e) 

i=0 

m«)e(!)(->)''e1 



i(Wt-l) 
i(iV t -l) 



(21) 



(22) 



(23) 



i=0 x y 1=1 

where substituting j v {y) in (20 ) leads to (21 1. Also, (22) is obtained by solving for In^u^^du in 
([21). Finally, ([23]) is given in [37]. 

Appendix II 
Proof of Theorem[2] 



Dropping the index k for the sake of convenience, the first term of ([T6j), T fc ( d ) is rewritten in terms 
of the channel directions, h as 



E ^ log 2 



|h r f full [ 2 

| h Tf|2 



E {log 2 (|h T f full | 2 )} - E {log 2 (|h T f| 2 )} . (24) 



Now, since |h T ff u n| 6 [0, 1], and |h T f | € [0, 1], we can safely assume that the angles Z(h c , f), Z(h c , f) 6 
[0,7r/2]. Using the triangle inequality for angles [38], we write 



Z(h c ,f) < Z(h c ,h c ) + Z(h c ,f), 
= Z(h,h) + Z(h c ,f), 



(25) 



where Z(h, h) = Z(h c ,h c ) is used to obtain (25). We can lower bound |h f| = cos (Z(h c ,f)) by 



cos 2 (Z(h c , f)) > cos(Z(h, h)) cos(Z(h c , f)) - sin(Z(h, h)) sin(Z(h c , f)) 
cos(Z(h,h))cos(Z(h c ,f)) 12 
cos 2 (Z(h,h))|h T f| 2 



(26) 
(27) 
(28) 
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where (27 1 is obtained by understanding that as the number of quantization bits increase, Z(h,h) will 
reduce and hence, sin(Z(h, h)) will be very small. The product sin(Z(h, h)) sin(Z(h, f )) can, therefore, 
be ignored. 



Substituting (28 1 in (24), we get 



E ^ log 2 



|h T f full | 2 

| h Tf|2 



< 



E{log 2 (|h T f full | 2 )}-E{log 2 (co S 2 (Z(h,h))|h r f| 2 )} 
-E{log 2 (cos 2 (Z(h,h))} (29) 



lo g2(e)^ 



i=0 



2-««,,(d 



j(JVi-l) 



(30) 



where ( |29j > was obtained using the fact that E {log 2 (|h T f full | 2 )} = E {log 2 (|h T f | 2 )} and (pOl), from 
Lemma Q] 

Appendix III 
Proof of Theorem [3] 

Dropping the index A; (and (k — 1)) for the sake of convenience, the second term of ( fT6] ), T^m, is 
upper bounded by 



E { log 2 



/Q(d)Q=|g f| +l 
P(d)«|g T ffuii| 2 + 1 



< E{log 2 (p (d) a|g T f| 2 + l) 



< log^l + z^aE^ff}) 

= log.Cl + ^aElUgflEllg^l 2 }) 

= ^(l + ^aiVtEllg^l 2 }) . 



(31) 
(32) 
(33) 



Here, ( pi) is obtained from Jensen's inequality. In (32 1 and (33 1, we use the relations g = ||g|| 2 g and 
^{||g|| 2 } = Nt, respectively. 

Now, the triangle inequality for angles is rewritten as [38], 

Z(g c ,f) > |Z(g c ,g c )-^(g c ,f)|, 
= |Z(g,g)-Z(g c ,f)|, 



since the angles are all positive and lie between [0,7r/2]. 
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Using trigonometry, sin 2 (Z(g c , f)) is lower-bounded by 

sin 2 (Z(g c , f)) > [sin(Z(g, g)) cos(Z(g c , f)) - cos(Z(g, g)) sin(Z(g c , f)) 

- 1 2 

-cos(Z(g,g))sin(Z(g c ,f)) 
cos 2 (Z(g,g)) 



(34) 
(35) 
(36) 



where ( |34| is obtained from the relations sin 2 (|0|) = sin 2 (0) and sin(i4 - B) = sm(A) cos(B) - 
sin(B) cos(A). Also, when sufficient quantization bits are used, sin(Z(g,g)) is extremely small. Hence, 
the product sin(Z(g, g)) cos(Z(g c , f)) can be ignored when the number of bits used for quantization is 
sufficiently large. Note that as cos(Z(g c ,f)) is close to (due to the generalized eigenvector relation 
that minimizes the value of |g T f|), sin 2 (Z(g c , f)) is very close to 1. This leads to ((36]) . 
We know that |g T f| 2 = cos 2 (Z(g c , f)). Using ([36]), we have 



E{|g J f| 2 } < E{l-cos 2 (Z(g,g))} 

N t 



where (37 1 is obtained from [36]. 



2 B f3 2 B , 



N t - 1 



(37) 



Appendix IV 
Proof of Theorem [4] 



We denote (from ( [19) ) 



A k = log 2 h + 2 Pfc , (i) 2Btnt _ fl> ( j + 2- B ^ log 2 (e) . 



(38) 



Note that A k is continuous and differentiable in B k ^ d y The partial derivative of A k in terms of -Bfc,(d) 
is obtained as 



dA k 2 Btat ~ Bk -<- d > 

2 Pk,(i) ( 2 B tot -B k ,, 



-Bk, 



dB, 



fc,(d) 



+ 2p Ul) + l) (l + 2 B — s,, ( d)) 



(39) 



Now, Afc is convex on the convex set C iff its gradient dA k /dB k r d \ is monotone, i.e. 

( ""^TiT " -^75T)1 [B h (d) (1) " B k (d) ^\ > . 
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From (38l, it is shown using simple algebra that if B^m) < B k ^ d ^ 2 \ then 





a) 



) 



< 



(40) 



2Stot-B fc , ( d) (2 » 



2 B t o t -B fc , (d) < 2 ) + 2 Pk>& + 1) (i + 2 



Btot — B, 




implying that the condition for convexity stated in this section is satisfied and A k is convex in B k ( d y 
Hence, the value of 5fc,(d) that minimizes will be a global optimal value and is obtained by equating 



([39]) to zero. A closed-form expression for the optimal Bwa) * s given by 



B k ,(d) = Btat ~ log 2 (l + p k ,(i) + ^ Pk,(i) 2Bt ° t+1 + (Pk,(i)) 2 ) ■ (41) 
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Fig. 1. Pictorial depiction of the Wyner model (described in Section ID). The solid line represents the desired signal, while the 
dashed line represents the interfering signal. 
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Fig. 2. The limited feedback model, described in Section [IV] to feedback quantized CSI of the desired and interfering channels 
using two separate codebooks. 



26 




6 8 10 12 14 

Received SNR of desired signal (dB) 



Fig. 3. Comparison of the actual and high SINR approximated sum-rates (in {4) and jTJ, respectively) for a — {0.001, 0.1, 1}, 
when N t = 4, p (d) = [0, 20] dB and K = 4 for full CSI. 
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Fig. 4. Sum-rates as a function of the number of cells (K) for a — {0.001,0.1, 1}, when Nt = 4, pM) — 10 d-B, for GEBF 
with full CSI. 
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Fig. 5. Mean loss in sum-rate for a = {10.001, 0.1, 1}, N t = K = 2, p (d) = 10 dB and B tot = 15. 




Fig. 6. Comparison of the mean sum-rates obtained using the proposed beamforming strategy, non-cooperative eigenbeam- 
forming and (cooperative) zero forcing, for p( d ) = 10 dB, N t = K — 2 and B to t = 6. 
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Fig. 9. Average data rate in a cell as a function of the received desired signal power for N t 



= K = 2 and B tot = 6. 



